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Abstract 

Parallel to the breakthroughs on the usage of simultaneous 
EEG-fMRI in neurocognitive studies and research, the 
occurrence of artefacts in the EEG signal induced within the 
fMRI scanner constitutes one of the challenges to be 
overcome in order to broaden the range of applications of 
such a technique. It is the case of the gradient artefact, 
provoked by the variation of gradient magnetic fields. 
Thereby, although a number of computational 
methodologies have yielded a satisfactory correction of the 
EEG signal, novel approaches have been proposed to 
improve the quality of the restored EEG. This work presents 
a novel proposal for modelling the variability of the gradient 
artefact template estimated during application of the 
established average artefact subtraction (AAS) method. 
Implementation of our proposed model and its combination 
with the AAS method allow an effective elimination of the 
gradient artefact from EEG recordings. Moreover, our 
approach performs an adaptive removal of the underlying 
artefact residuals, according to the signal slope parameter 
characteristics. 
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Introduction 

The usage of simultaneous EEG-fMRI holds a large 
number of applications in recent years because its 
powerful capability of providing new insights into the 
brain function (Villringer et al., 2010). Albeit such a 
technique was firstly applied in the field of epilepsy 
(Warach et al., 1996; Seeck et al., 1998), over the years 
several other neuroscientific studies, research, and 
applications which make use of the combined EEG- 
fMRI have been proposed and reported in the 
literature (Villringer et al., 2010; Ritter and Villringer, 
2006). 

However, enhancements still remain to be accomplished. 


regarding the induction of artefacts in the EEG signal 
by the magnetic fields of the fMRI scanner. In this 
respect, the quality of the EEG signal acquired during 
continuous fMRI depends on the performance of 
hardware solutions as well as computational 
approaches to minimize and remove the common 
types of artefacts arising in the electroencephalogram 
(Allen et al., 1998). That is the case of the so-called 
"gradient" or "imaging acquisition" artefact (Allen et 
al., 2000; Mulert and Hegerl, 2009; Ritter et al., 2010). 

Gradient artefacts occur in the EEG signal induced by 
the application of rapidly varying gradient magnetic 
fields for spatial encoding of the magnetic resonance 
signal (MR), and radiofrequency pulses (RF) for spin 
excitation in the circuit formed by the electrodes, leads, 
patient, and amplifier. They have a waveform which is 
approximately the differential waveform of the 
corresponding gradient pulse (Anami et al., 2003; 
Ritter et al., 2010). The amplitudes of imaging artefacts 
can be several orders higher than the magnitude of the 
neuronal EEG signal. According to Ritter et al. (2010), 
artefacts associated with the gradient switching have 
amplitude at 10 3 to 10 4 pV, whereas those arising from 
RF have magnitude up to around 10 2 pV. In turn, the 
bandwidth of the standard clinical EEG is exceeded by 
the gradient artefact frequencies, in addition to being 
overlapped by discrete harmonic artefact frequency 
intervals, or "frequency bins". As discussed by Niazy 
et al. (2005), the fundamental of each respective 
frequency bin corresponds to multiples of the inverse 
of the EPI slice time parameter. 

In order to minimize the induction of gradient 
artefacts in the EEG signal, a number of hardware 
solutions have been suggested (Anami et al., 2003; 
Ritter et al., 2010; Mullinger et al., 2011). Concerning 
post-processing artefact correction methods, frequency 
and time-domain techniques have been proposed as 


55 



www.seipub.org/spr 

well, like frequency-domain filtering (Hoffmann et al., 
2000; Sijbers et al., 1999; Sijbers et al., 2000), average 
artefact subtraction (AAS) (Allen et al. 2000), principal 
component analysis (Negishi et al., 2004; Niazy et al., 
2005), independent component analysis (Mantini et al., 
2007), and spatial filtering (Brookes et al., 2008). 

The average artefact subtraction (AAS) is the most 
established computational methodology for imaging 
artefact correction. It consists of the calculation of an 
average artefact waveform (or template) over a fixed 
number of samples, which is then subtracted from the 
EEG signal at each sample (Allen et al., 2000). The 
usage of the AAS method allows a substantial 
attenuation and removal of the gradient artefact. 
Nevertheless, as it is assumed that the artefact 
waveform possesses a constant morphology over the 
time, and does not take into account its variability, 
artefact residuals arise in the EEG signal after the 
template subtraction. Further, effectiveness of the AAS 
depends on accurate sampling of the gradient artefact 
waveform and timing- alignment of the fMRI 
equipment and EEG clocks (Mandelkow et al., 2006; 
Goncal ves et al., 2007; Freyer et al., 2009). Thereby, 
additional restoration methods must be carried out in 
order to achieve a satisfactory signal correction (Allen 
et al., 2000; Niazy et al., 2005, Gongalves et al., 2007). 

The current work presents a novel approach for 
modelling and removal of the gradient artefact 
residuals resulting from application of the AAS 
method. In this sense, a specific model is proposed to 
evaluate and quantify the variability of the artefact 
template waveform. Our approach is modified from 
Ferreira et al. (2012), and based upon the information 
of variability contained in the difference between 
consecutive samples of the digital EEG signal 
(Cluitmans et al., 1993; Van de Velde et al., 1998; 
Ferreira et. al, 2013a). 

Materials and Methods 

Patients 

The EEG recordings were collected simultaneously 
with fMRI data for a research focused on epilepsy and 
post- traumatic stress disorder (PTSD) (Van Liempt et 
al., 2011), jointly developed by the department of 
Psychiatry of Universiteit Medisch Centrum Utrecht, 
the Research Centre Military Mental Health Care in 
the Dutch Central Military Hospital in Utrecht, and 
the Department of Research and Development of the 
Kempenhaeghe Epilepsy Center in Heeze, The 
Netherlands. 
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The data were recorded from military veterans with 
PTSD which were in mission abroad through the 
outpatient clinic of the Military Mental Health Care. 
All participants were male and aged between 18 and 
60 years. EEG signal recordings from 15 subjects were 
used in order to apply and evaluate the methodology 
we have proposed in this work. 

Protocol for Data Acquisition 

A 3 T Scanner (Philips, Eindhoven, Netherlands) was 
employed to carry out the functional magnetic 
resonance imaging scanning. An MRI-compatible 64 
channel polysomnograph was used to collect one ECG 
channel, two EOG channels, one EMG and 60 EEG 
channels. The EEG electrodes positioning was in 
accordance with the 10 - 20 international system 
electrodes placement. The sampling frequency for EEG 
signal acquiring was 2048 Hz. 

Simultaneous EEG-fMRI data were collected during 45 
minutes, after application of the EEG cap. The subjects 
were scanned using a functional echo-planar imaging 
sequence with 33 transversal slices (thickness 3 mm, 
TE 30 ms, and TR 2500 ms). 

General Overview of the Average Artefact Subtraction 
Methodology 

The basic idea of the average artefact subtraction 
method consists of estimating an average template of 
the gradient artefact waveform observed in the raw 
EEG signal. Such a template is then subtracted from 
the electroencephalogram at those regions where the 
artefact occurs (Allen et al., 2000). The average artefact 
subtraction can be described by Eq. (1) (Ferreira et al., 
2012 ): 

EEG corcti = EEG 

raw,i F i s (i) 

where EEGcorct and EEGraw are, respectively, the 
corrected EEG after template subtraction and the raw 
EEG signal; T corresponds to the template to be 
subtracted; i runs over the number of samples within 
EEG signal; and is rims over the template segment. 

According to Allen et al. (2000), in order to calculate 
the template T, the EEGraw excerpt is divided into 
segments of length S. The average of samples situated 
at the same template position is corresponds to the 
template sample 73 . For continuous fMRI acquisition, 

S is defined as the EPI slice time parameter (ST). 
Because of the misalignment between the EEG and 
fMRI clocks, interpolation and extrapolation within 
the segments must be performed as well (Allen et al., 
2000 ). 
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Overview of the Average Artefact Subtraction 
Methodology Employed in This Work 

In this work, the AAS method was also carried out by 
using the model of Eq. (1). For estimation of the slice 
time (ST), we used the methodology proposed by 
Garreffa et al. (2003). Observation of the raw data 
recorded within the MR scanning revealed the 
occurrence of the typical peaks corresponding to the 
onset of each MR slice, in such a way that ST could be 
estimated by evaluation of the time interval between 
two subsequent peaks. For the data under analysis, the 
value of ST was calculated at 155 ± 1 samples, which 
corresponded to the time interval of 75.68 ± 0.50 ms 
(Ferreira et al., 2012). 

Hence, the raw EEG was divided into segments of 155 
samples for estimation of the artefact template. In case 
of segments which contained more or less than 155 
samples, extrapolation or interpolation was done in 
order to compensate the misalignment between the 
EEG and fMRI clocks. The number of segments 
considered for average was 32 slices. 

To implement the average artefact subtraction method, 
we used the approach described in Ferreira et al. 
(2012). 

Modelling of the artefact waveform variability 

To identify and remove the remaining residuals 
resulting from application of the AAS method, a 
model for quantification of the artefact waveform 
variability has been proposed, modified from Ferreira 
et al. (2012). As discussed by these authors, estimation 
of the template variability could be adaptively 

implemented by multiplying each artefact template 

sample by an estimated factor di, and then the 
subtraction indicated in (1) is performed as follows: 

EEG corcti = EEG 

raw,i Wi, ( 2 ) 

Eq. (2) is, therefore, an adaptive extension of Eq. (1), 
whose limit referred to as the filter parameter Si allows 
refining the value of EEGcorct (Klados et al., 2009; 
Ferreira et al., 2012): 

lim EEG corcti = true EEG, (3) 

dj -Adj 

in such a way that, as the value of Si is adapted and 
approximates to the optimal parameter value at, the 
value of EEGcorct, i would tend to the true EEGi. 

In the current work, a was changed to a new 
parameter a', and Eqs. (2) and (3) have been further 
modified, resulting in: 
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E E G C orct,i E E G corc t i CL iRi 

, r 

EEG corct i+1 E EG C0rC f :i i+i + cl iRi (4) 
where the product of the corresponding elements a' 
and R acts as an initial estimation of the remaining 
residual in the EEGcorct. Thus, Eq. (4) allows calculating 
a refinement for the EEG signal corrected by the AAS 
method. And the respective restored EEG (EEGcorct) is 
calculated by subtracting the estimation of the 
remaining residual, S'iRi, from the corresponding 
unrefined values within EEGcorct. 

As the product of the corresponding elements of a' 
and R should model the variability of the gradient 
artefact, for their estimation we have sought and taken 
into account specific properties of the raw EEG, the 
EEG corrected by AAS, and the true EEG signal which 
could reflect such a variability. In this way, the signal 
slope parameter was assumed to possess such 
characteristics. 

The analysis of the signal slope represented by the 
difference between consecutive samples can be 
employed to detect high-frequency properties and 
artefacts in digital signals (Barlow, 1983; Cluitmans et 
al., 1993; Van de Velde et al., 1998; Ferreira et al., 2012). 
The use of the signal slope is in accordance with 
Cluitmans et al. (1993) and Van de Velde et al. (1998) 
who described such a signal parameter associated 
with the large signal magnitudes as being particularly 
useful for detection of the high-frequency properties of 
the muscle artefact. In the same way, observation of 
the EEG recorded within the fMRI scanner reveals that 
large slopes and magnitudes of the signal can be 
attributed to the gradient artefact and its residuals as 
well (Koskinen and Vartiainen, 2009; Ferreira et al., 
2012). Thus, the slope variation was also employed to 
identify the imaging acquisition artefact, described as 
follows. 

According to Ferreira et al. (2012), the larger slopes 
related to the sharp wave activity of the gradient 
artefact residuals are used to identify whether the EEG 
samples are artefact-free or not. In order to perform 
this identification, a slope threshold (tins) is estimated 
in such a way that if the sample has signal slope larger 
than thrs, it is then identified containing artefact 
interference, thrs can be estimated, for example, 
considering the probability distribution of the signal 
slope (Cluitmans et al., 1993; Van de Velde et al., 1998). 

In this way, we have taken into account the probability 
distributions of an artefact free EEG (EEGtme) and the 
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EEGcorct excerpts, both ones picked up from the same 
EEG channel of one specific subject. The true EEG was 
obtained from a raw EEG interval recorded within the 
scanner with no gradient magnetic field variation. It 
was observed that more than 90% of the values for the 
signal slope of the EEGcorct, diff (EEGcorct), are above the 
threshold calculated using the EEGtme excerpt (Ferreira 
et al., 2013b): 

thrs = /t d j //(EEGtrue ) + 3 cr di//(EEGtrue ) (5) 

Furthermore, the probability distribution of the slope 
parameter estimated for the true EEG (EEGtme) closely 
resembles Gaussian distributions (Cluitmans et al., 
1993), in such a way that the confidence interval 
associated with diff <EEc,„ e ) + diff(EEG„„> encompasses 
approximately 99.5% of the distribution (Papoulis and 
Pillai, 2002). Thus, estimation of thrs was performed by 
using Eq. (5). 

The maximum absolute value of the difference 

between consecutive samples of the signal EEGcorct 
corresponds to the parameter n which is related to 
diff (EEGcorct) by the following expression (Ferreira et 
al., 2013a): 

r t = max \diff (EEG corct )\ (6) 

where i is the subscript of the maximum slope within 
diff (EEGcorct). 

The two consecutive samples EEGcorct, i and EEGcorct, i+i 
associated with n are adapted as follows: 

EEG corcti £ , £ , fj C0rC ( j ; Lj 

EEG corct i+ i EEG corct i+1 -b L; (7) 

where 

Li —p Xr t (8) 

p is the adaption factor (0 < p < 1) to be applied to the 
parameter n. We experimentally determined and set 
the value of p as 0.925 (Ferreira et al. 2013b). Therefore, 
the product ai'Ri of Eq. (4) corresponds to Li in Eq. (8). 
In Eq. (7), the sign of Li is set positive when 
EEG cord, i ^ EEGcorct, i+lf and vice-versa. The signal EEGcorct 
in equation (6) is then replaced by the modified signal 
EEGcorct, which contains the adapted samples EEGcorct, ; 
and EEGcorct, i+i. Eqs. (6), (7), and (8) are iteratively 
recalculated until n < thrs (Ferreira et al. 2013b). 

According to our approach, therefore, in addition to 
acting as an indicator of the presence of the gradient 
artefact residuals in the EEGcorct (Cluitmans et al., 1993; 
Van de Velde et al., 1998), the signal slope parameter, 
diff (EEGcorct), is also used to estimate the amount of 
correctness (within Li) associated with the template 
variability, to be applied and restore the signal EEGcorct, ; 
(Ferreira et al., 2012). Thereby, removal of the artefact 
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is performed in an adaptive way by approximating the 
value of the parameter n to the slope characteristics of 
the raw EEG, corresponding to the threshold thrs. 

In comparison to the method presented in Ferreira et 
al. (2012), the usage of Eqs. (5) - (8) shows itself to be 
more robust and effective on removing stronger 
gradient artefact interference. Although both 
approaches make use of the information of variability 
contained in the signal slope parameter, the former is 
based upon normalized values of the signal slope 
parameter, whereas according to the latter, the 
magnitude of the difference between consecutive 
samples is used instead (Eq. 6). Moreover, in this work 
we did not take into account the standard deviation 
associated to the template average for modelling the 
variability of the residuals. It was noticed that 
estimation of such variability is more accurate as 
performed by Eqs. (5) - (8). 

Validation 

A set of validation EEG channels, FP1, F3, F8, T5, P3, 
Oz, CPI, FC5, AFz, F6, C2, TP7, CP4, and POz, was 
used for evaluation of the proposed model for the 
artefact residuals from the AAS method and respective 
EEG restoration. In this sense, the median attenuation 
of the power spectrum associated with the artefact 
frequency bins (Niazy et al., 2005) as well as within the 
range of the clinical EEG (0 - 24 Hz) was assessed. 
Comparison with the signal restoration achieved by the 
usage of a linear low-pass filter was also performed. 

Results 


Raw EEG 



FIG. 1 REPRESENTATIVE RAW EEG SELECTED FOR 
ILUSTRATION OF APPLICATION OF OUR PROPOSED 
METHODOLOGY 

Fig. 1 depicts a typical EEG excerpt containing strong 
artefact interference (at 3.7 mV pk-pk), selected from 
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the recordings of one subject, electrode position F6. 
The steep wave activity associated with the gradient 
artefact can be observed in the EEG signal, which is 
completely obscured by the gradient artefact. Such an 
EEG excerpt was segmented as well as interpolated or 
extrapolated depending on the need, according to the 
AAS method (Allen et al. 2000). The set of 32 slices 
used for template averaging is shown in Fig. 2. 
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FIG. 2 SET OF 32 SUBSEQUENT SEGMENTS (32 SLICES) IN 
WHICH THE SIGNAL OF FIG. 1 WAS DIVIDED FOR ARTEFACT 
TEMPLATE ESTIMATION (155 SAMPLES) 

Corrected EEG by application of the AAS method 

1500 1 1 1 1 1 


1000 - 
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FIG. 3 CORRECTED EEG SIGNAL (EEGcorct), RESULTING FROM 
THE APPLICATION OF THE AAS METHOD 

Restored EEG (cleaned from residuals) 



FIG. 4 RESTORED EEG SIGNAL (EEGcorct), RESULTING FROM 
THE APPLICATION OF EQS. (5) - (8) IN THE SIGNAL EEGcorct 
DEPICTED IN FIG. 3 
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Fig. 3 shows the signal EEGcorct resulting from the 
application of the average artefact subtraction 
methodology in the raw EEG signal depicted in Fig. 1. 
A considerable amount of residual artefact remained 
in the EEGcorct after the template subtraction (at 2.0 mV 
pk-pk). It can be attributed to the large amplitudes and 
variability of the artefact waveform over the time, as 
observed in Fig. 2 for the segmented EEG slices. Hence, 
in order to restore the EEG signal, the model described 
by Eqs. (5) - (8) was then applied to the signal EEGcorct 
of Fig. 3. The resulting EEGcorct is depicted in Fig. 4. 
From this figure, it can be observed that the artefact 
residuals were strongly cleaned up. 

As shown in Figs. 5 and 6, our approach removes 
high-frequency artefact residuals components, as it is 
performed by a low-pass filter. A linear low-pass filter 
was used for comparison purposes, set up as a 55- 
coefficients FIR, cut-off frequency 70 Hz. Moreover, 
our approach also causes an adaptive attenuation of 
underlying artefact residuals components within the 
EEG bandwidth, which can be noticed in Fig. 6d. 

According to Niazy et al. (2005) underlying gradient 
artefact components arise and leak into the frequency 
bins intervals, region of ±1 Hz around the multiples of 
the inverse of ST. The linear low-pass filter does not 
yield effective elimination of frequencies below the 
filter cut-off frequency and, therefore, it also does not 
attenuate the respective underlying artefact residuals 
in the EEG bandwidth, as depicted in Fig. 5b (short 
peaks approximately even spaced) and Fig. 6b. On the 
other hand, low-frequency artefact components are 
attenuated by application of our approach in such a 
way that those residuals have been cleaned up in the 
restored EEG signal of Fig 5a. 

In Table 1, the median spectral power attenuation 
provoked by application of our approach and the 
linear filter is shown. The attenuation was calculated 
using Eq. (9): 

Attenuation = 100 x absolute 5r ) % (9) 

where Saas is the power spectrum of the EEG corrected 
by the AAS method; and Sr is the power spectrum of 
the restored EEG. Regarding the EEG validation 
channels, the median spectral power attenuation 
within the frequency bins 13.2, 26.4, 39.6, 52.9, 66.1, 
79.3, and 92.5 Hz achieved by the application of our 
approach is higher than that provoked by the linear 
filter for frequencies below the filter cut-off frequency. 
It even can be noticed that application of the filter 
amplifies the energy of some frequency bins. Our 
approach has caused higher power attenuation in the 


Set of 32 slices used for template estimation 
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ranges 4-8 (6.6%) and 8-12 Hz (17%) as well. The EEG after application of Eqs. (5) - (8) was 2.9 pV pk- 
median value of the remaining artefact in the restored pk. 


a Restored EEG by the model described in Eqs. (5) - (8) 



Time (s) 


b Restored EEG by linear filtering 



FIG. 5 EEG RESTORATION ASSOCIATED WITH THE SIGNAL EEGorct DEPICTED IN FIG. 3, OBTAINED BY (A) THE PROPOSED MODEL 
FOR ARTEFACT RESIDUALS DESCRIBED IN EQS. (5) - (8), AND BY (B) A LINEAR LOW-PASS FILTER 






FIG. 6 COMPARISON BETWEEN THE POWER SPECTRUM OF THE EEG CORRECTED BY THE AAS METHOD (FIG. 3) AND: (A) THE 
POWER SPECTRUM OF THE EEG RESTORED BY OUR APPROACH (FIG. 5A); (B) THE POWER SPECTRUM OF THE EEG RESTORED BY 
LINEAR FILTERING (FIG. 5B); AND (C) THE POWER SPECTRUM OF THE RAW EEG SIGNAL (FIG. 1). (D) POWER SPECTRUM OF THE 
RESIDUAL NOISE RESULTING FROM APPLICATION OF OUR APROACH AND LINEAR FILTERING IN THE SIGNAL OF FIG. 3 
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TABLE 1 MEDIAN SPECTRAL POWER ATTENUATION PROVOKED BY APPLICATION OF OUR APPROACH AND LINEAR FILTERING AFTER AVERAGE TEMPLATE 

SUBTRACTION, CONSIDERING THE EEG VALIDATION CHANNELS 


Restoration 

Method 

Attenuation (%) j 

EEG Frequency Bands (Hz) 

Artefact Frequency Bins (Hz) j 

0-4 

4-8 

8-12 

12-24 

13.2 

26.4 

39.6 

52.9 

66.1 

79.3 

92.5 

Our approach 

1,1 

6,6 

17,0 

34,7 

28,6 

56,4 

72,9 

83,6 

87,9 

94,8 

96,1 

Linear filtering 

2,7 

4,9 

6,7 

37,7 

26,5 

184,2 

297,3 

228,6 

264,2 

98,2 

93,4 


Discussion 

Application of the average artefact subtraction method 
allows substantial removal of gradient artefact 
components. However, partial failure of such a 
method occurs because the variability of the gradient 
artefact waveform over the time, as can be observed in 
Fig. 3. Thereby, the usage of further correction 
approaches is required for removal of the artefact 
residuals (Allen et al., 2000; Niazy et al., 2005). 

High-frequency artefact residuals which exceed the 
range of the clinical EEG can be promptly cleaned up 
by applying linear low-pass filters, as depicted in Figs. 
5b and 6b. However, the underlying residuals within 
the EEG bandwidth are not removed by linear filtering, 
demanding further correction techniques (Allen et al., 
2000). In this respect, we have proposed the model of 
Eqs. (5) - (8) to quantify and remove the remaining 
artefact residuals. Such a model constitutes a non- 
linear filter approach based on the signal slope 
parameter, which attempts to minimize the variability 
over the time associated with the estimated artefact 
template. 

In addition to eliminating the high-frequency artefact 
residuals components, our approach performs an 
adaptive attenuation of the underlying artefact 
residuals, as depicted in Figs. 4, 5, and 6, and observed 
in Table 1. Unlike a linear low-pass filter, our 
proposed model adapts only the samples of the signal 
which have signal slope higher than thrs (Ferreira et al., 
2013a). Hence, the other samples are not adapted and, 
in consequence, are not distorted by the filtering 
process. This fact explains the lower signal distortion 
provoked with the application of our approach. 
Further, as the model described in Eqs. (5) - (8) takes 
into account a specific property of the artefact (i.e., the 
signal slope parameter), it attenuates the residuals in 
accordance with the artefact slope characteristics. 
Thereby, application of Eqs. (5) - (8) in the signal 
EEGcorct adapts firstly the signal slope parameter 
associated with the underlying artefact residuals, 
which are larger than that of the clinical EEG. 


Therefore, artefact residuals and eventual leakage 
around the corresponding frequency bins (Niazy et al., 
2005) are attenuated without a preliminary need for 
using additional correction techniques. 

It is noteworthy that our methodology is even effective 
within a strong artefact interference scenario, as 
depicted in Fig. 1, whereas the method presented in 
Ferreira et al. (2012) is suitable to remove artefact 
residuals with lower order magnitudes (at 150 pV pk- 
pk). To smooth the restored EEG signal, a moving- 
average filter could be used as well after application of 
Eqs. (5) - (8) (Ferreira et al., 2013b). In future work, 
further improvements concerning the computational 
performance of such a non-linear filter shall be done. 
Also, our methodology shall be assessed for gradient 
artefacts induced in the EEG signal within other types 
of fMRI scanners. Finally, the proposed approach shall 
be evaluated for digital filtering and correction of 
other types of artefact as well. 

Conclusions 

In this work, a non-linear filter approach is proposed 
for modelling of the gradient artefact residuals 
resulting from application of the established AAS 
method. 

Implementation of our approach is based upon the 
idea of minimizing the variability associated with the 
artefact template waveform. In this respect, the signal 
slope parameter, defined as the difference between 
consecutive samples of the digital EEG signal, is 
employed to estimate such variability. 

The proposed model performs an adaptive artefact 
filtering of high-frequency artefact residuals 
components as well as underlying components in the 
EEG bandwidth. Therefore, combination of our 
modelling approach with the AAS method allows 
effective artefact residuals attenuation and EEG 
restoration. Hence, our method could be used and 
even combined with other artefact correction 
techniques in order to improve the quality of the EEG 
signal acquired within the fMRI scanner. 
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